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ABSTRACT 


A general theory of the diffusion cooling of neutrons in a finite solid moderator assembly 
is given. In the present theory, starting from the transport equation in the diffusion ap- 
proximation, the energy-dependent part of the solution is expressed in terms of generalized 
Laguerre polynomials of order unity and degree n, each term of which is weighted by a Max- 
wellian distribution function. In the approximation in which only the first two terms of the 
Laguerre expansion are retained, one obtains exactly the same expression for the decay con- 
stant of the fundamental mode as in the variational approach of Nelkin and of Singwi and 
Kothari. In the next higher approximation, one gets a correction which we call the “non- 
Maxwellian” correction to the diffusion cooling constant, which correction is nearly twenty to 
twenty-five per cent. The merit of the present theory is that it avoids the assumption of a 
neutron temperature which was the weak point of the variational approach. Expression for 
the correction to the diffusion cooling constant arising from terms of order B® has also 
been derived and it turns out that this correction has a sign opposite to that of the non- 
Maxwellian correction and in magnitude it is quite significant, depending on the value of the 
buckling, B? of the assembly. Non-diffusion correction to the decay constant has also been 
estimated and it is found that it is negligible compared to the multi-velocity correction. Ex- 
pressions for the mean energy of neutrons inside and outside the assembly are derived. Finally 
the variation of the diffusion cooling constant with the temperature of the moderating as- 
sembly is studied in the two limiting cases: (i) 1/@—> co and (ii) 6/7’ > oo, 6 being the Debye 
temperature of the solid. Exact evaluation of this variation is shown graphically in the case 
of beryllium. Comparison of the theoretical values of the diffusion cooling constant with the 
experimental values for beryllium, beryllium oxide and graphite has been made. If the trans- 
port mean free path is assumed to be independent of energy, the agreement between theory 
and experiment is good in the case of graphite and beryllium, but bad in the case of beryllium 


oxide for which the experimental value seems to be rather high. 


1 Based on work performed under the auspices of the U.S. Atomic Energy Commission. 
2 On leave of absence from the Atomic Energy Establishment, Trombay, Bombay. 
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In an earlier paper, hereafter referred to as I, Singwi and Kothari ’(19% 
had derived an expression (Eq. 8 of I) for the decay constant, / of the funda- 
mental spatial mode of the neutron flux, as a function of the geometrical buck- 
ling, B® of the moderating assembly. It is generally expressed in the form 


A= av + D, B’—CD, Br, 


where the various symbols have their usual meaning. This expression is valid — 


for times long compared to the slowing-down time of neutrons, which are in- 
jected into the assembly in the form of a pulse. The significance of the first 
two terms in the expression for 4 is very well understood and there exists a 
fairly good agreement between theory and experiment as regards the value of 
the diffusion coefficient D,. This, however, is not the case with the last term 
in which the coefficient C is called the diffusion cooling constant. In a finite 
assembly there is a preferential leakage of fast neutrons; the rate of leakage — 
depends on the size of the assembly and on the energy dependence of the 
diffusion coefficient. As a consequence of this, the neutrons inside the as- 
sembly cool below the moderator temperature. Thermal equilibrium will, how- 
ever, be established as a result of the transfer of energy between the neutrons 
and the solid. The rate at which this equilibrium is established will depend on 
the efficiency of the mechanism of energy transfer. In fact the constant C is 
a measure of this efficiency. At low temperatures of the moderating assembly, 
the energy transfer is so slow that a thermal equilibrium may never be at- 
tained. 

deSaussure and Silver (1958) have tabulated the experimental values of C 
as determined by different workers for beryllium at room temperature. It will 
be seen from their table that the values of C range from 0 to 3.05 +0.65 em?. 
Their own value is 1.1+0.9 em?. The experimental position! is, therefore, very 
unsatisfactory. The theoretical value of C as calculated by Nelkin (1958) is 
apparently in good agreement with the experimental value of deSaussure and 
Silver. However, in a more realistic theory if one takes into account the varia- 
tion of the transport mean free path with neutron energy this agreement is 
spoiled (paper I). 

The earlier theories of Nelkin and of paper I were based on the assumption 
that the departure of the neutron energy spectrum from the equilibrium Max- — 
wellian spectrum was small; and that one could associate a temperature 7’, 
(if it exists!) with the neutrons. A trial function of the Maxwellian form for 
the neutron spectrum was taken with the neutron temperature 7’, as a varia- 
tional parameter. The Rayleigh-Ritz variational principle was used to determine 
the lowest eigenvalue of the decay constant. Nelkin examined the accuracy of 
the assumption of a neutron temperature by comparing his calculation with the 
exact results of Hurwitz et al. (1956) for the case of a heavy gaseous moderator. 


a For a general experimental discussion of the diffusion cooling constant see the review 
article (Progress in Nuclear Energy, Series 1, Physics and Mathematics, Vol. 2. Pergamon 
Press, London, 1958) by von Dardel and Sjostrand. 
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3 His estimate was that the error in the diffusion cooling constant was of the 
_ order of twenty-five per cent. If the moderating assembly is small, there will 
be a large leakage of neutrons and consequently it becomes doubtful whether 


_ the association of a temperature to the neutrons inside the assembly has any 


~ pletely avoided. Mention may be made here of a numerical calculation, based 


_ degree n, each one of which is weighted by a Maxwellian distribution function. 
_ Such an expansion has also been used recently by Kazarnovsky et al. (1958). 


physical meaning at all. The main purpose of this paper is to develop a more 


_ .general theory in which any such assumption of a neutron temperature is com- 


on a heavy gas model by Beckurts (1957), who found that considerable devia- 
tions from a Maxwell spectrum occur even for rather small bucklings. 

The starting point of the theory is to expand the velocity distribution of the 
neutron flux in terms of generalized Laguerre polynomials of order unity and 


The number of terms required in the expansion will obviously depend on the 
closeness of the actual spectrum to the Maxwellian spectrum. On retaining the 


first two terms of the Laguerre expansion, we get the same expression for C 


_. as obtained in I using the variational principle. The third term in the expan- 


sion gives a correction (to the order 5*) to the diffusion cooling constant which 
is estimated to be nearly twenty per cent. We would call this as the ‘“non- 
Maxwellian” correction. 

For an accurate experimental determination of C it is necessary that the 
moderating assembly be small, i.e. its buckling, B? be large. In such a case it 
is not possible to neglect terms of order B® in the expression for the decay 
constant. We have derived expressions for terms of order B® and B® in the 
approximation in which only the first two terms in the Laguerre expansion are 
retained. Expression for the term of order B® has also been obtained when the 
first three terms in the expansion are retained. It turns out that this correction 
has a sign which is opposite to that of the non-Maxwellian correction and_ it 
increases in magnitude at low temperatures. For the smallest beryllium assembly 
used by deSaussure and Silver, it is nearly fifty per cent of the non-Max- 
wellian correction to the diffusion cooling constant. 

The magnitude of the “non-diffusion” correction to the cooling constant has 
also been estimated and it is found to be negligible compared to the non- 
Maxwellian correction. 

Expressions for the mean energy per neutron in the flux inside and outside 
the assembly are given, which clearly indicate the effect of the variation of the 
transport mean free path with neutron energy. The inside energy, in general, 
differs from the outside energy; it is the latter which is measured experiment- 
ally. 

The temperature dependence of the diffusion cooling conssant has been stud- 
ied in the two limiting cases: (i) 7’/9—0o and (ii) T/6 +0, 0 being the 
Debye temperature. In the latter case C varies as (0/T)’, which indicates that 
for very low temperatures the diffusion cooling constant will be so large that 
thermal equilibrium may not be possible at all, whereas in the former C is 
constant. A graph showing the calculated variation of C with (0/ T) is also 
given. An experimental verification of this result would be interesting. 
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2. Mathematical formulation 


The Boltzmann equation for the angular flux f(r, #, Q, t) is 


Eogith 4,8) ~Q.- grad f(r, E, Q, t)—f (r, E, Q, t)[S.(t, B, Q)+ Sa(t, B)) 
Vv 


+9(t, 8, 9)+ [ dE"dQ! I, B, 2,0 3.(0, BE, Q'>Q) (1) 


(see Weinberg and Wigner (1958), p.. 223), and where the symbols have their | 
well-known meaning. v is the neutron velocity corresponding to energy E. Since ~ 
our interest is in developing a theory for the pulsed neutron experiments, we 
shall take the source term S(r, #,Q) to be zero. In the P, (or the diffusion) 
approximation equation (1) takes the form 


1 aO(r,£,t) In(Z 


1524 (r, H, t)+3a(E) b(t, By t)+$(r, B, t) 5, (EB) 


v ot 3 c. 
- | 60, E’,t)3.(E'> E)dE’=0, (2) 
0 
where $(t, B, t)= | dQi(r, E,Q. 0). 


The validity of Eq. (2), as is well known, is based on two approximations: 
(i) the neglect of the second and higher terms in the angular distribution of 


f and (ii) the neglect of the al term in the expression for J, where J is the 


neutron current. We shall examine the influence of approximation (i) when we 
discuss the magnitude of the non-diffusion correction to the diffusion cooling 
constant, and as regards the second approximation, it is easy to show that it 
is valid so long as the scattering cross section is much greater than the ab- 
sorption cross section. 

By making a Laplace transform, the time variation can be separated from 
the transport Eq. (1), and the problem reduces to a time-independent equation, 
in which the time constant 4, or rather 4/v appears as an eigenvalue. This 
eigenvalue equation was solved in the diffusion approximation in paper I by 
using the variational method, assuming that the space and energy dependence 
of f can also be separated out. 

Our starting equation is (2). Let us introduce a dimensionless variable x= v/v; 
where v=(2k)7'/m)!, T’ is the. moderator temperature, ky the Boltzmann’s 
constant and m the neutron mass. Equation (2) then becomes 

oo 
oi) (tae. t) ae 
Va £ ) 3 le, (x) V* b (8, a, t) + Da (x) H(r, 2, t) oe th Unstse &) | >s(@—> 2’) dz’ 
0 


= { f(r, x’, t) S,(x’>2x)da’=0, (3) 


where in (3) we have put v equal to unity. 
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In Eq. (7) J denotes the transport mean free path and where we have omitted 


the suffix tr. 
Let us now expand n;,(r,t) in terms of the eigenfunctions of the equation 


V? Xi (r) + Bi Xi (r) =0, 
where B? is the buckling corresponding to the Ith harmonic. Now 
or (tr, N= 2 st (0) (er (10) 
Since we are interested in times sufficiently long compared with the slowing- 


down time, all higher harmonics must have died out and we shall be left with 
only the fundamental mode. We then have 

mg (r, 1) =8 (0) Lor), (11) 
where Xo (r) is the eigenfunction corresponding to the lowest eigenvalue B?. 
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for k=1 and i=0 and 1. 
Let the solutions of (12) be 


gO = e4t ) 
rt 


and s® = ge7 
; ‘ms 


Substituting (13) in (12) we have 


1 it 1 1 | 
2a (ees Seay (ze ay ay ote is P = 2 ae f 
(;2)* (7a) 43 fa +3 Yoo + 0 (S ina a=0 (14a) 


x 


and 


me ee) A- ial Aa + 5 ly B+ 5 ly Bat (Sahota (San tay =O. (14 b) 
00 11 


Eliminating « in the above two equations, we get the following quadratic equa: | 
tion in A 
4 4B WStt.e 8 
#+21|—K- = are Giga wb egsbay'n K 
3Vn°"  9Va\V2 8 ™ 30m" 
UG ai jn 16 16 B’ 
—~——= | eee ee peace —2)= 
yale, g “00 I) +958 loo Yur + 77 (Li loo — 41) = 0, (15) 


where we have made use of the following: 


which can be easily deduced using Eqs. (8) and (5). We also note that fire = fret» 
The macroscopic absorption cross section Sz is assumed to be inversely propor- 


tional to velocity 7.e., oe 
es 
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sl bl a)ne d= 7? j Air (EB) E exp (— E/k, T) a E, 


_ which is the same as A, of Eq. (9) of paper (I). 


Also a= 7 Al ler (a) @e-* LO? (a) LY (w) dx 


= 


“al bi, (2) we~* (2—ax) da 


a 1 

eS =V2A)—= Ay, (20) 

vor 

, where A= Goan | 1, (E) E'** exp (— E/ky)T) dE. (21) 

4 0 

5 It is shown in the appendix (Eq. A-20) that 

g 4y,= My, (22) 

_ where M, is given by Eq. (4) of paper I. Using the above results (18) becomes 

4 123,045? Dj— BD, 0 (23) 
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Vx Dy (2 -3) (24) 
where aay 64 re ) AS 


which is the same as equation (12) of paper I except that one of the D,’s 
here is taken with B*. Thus if we retain terms corresponding to 1=0,1 and 
k=0,1 in Eq. (7) we get the same result as was obtained earlier using the 
variational principle. To get a non-Maxwellian correction to the diffusion cooling 


OT 


constant CO, we must carry out the next approximation, 7.e., retain terms cor- — 


responding to +=0.1,2 and k=0,1 and 2. This will be done in Section 6. 


4. Terms of order B® and B® 


Since the moderating assemblies have to be necessarily smal] for an accurate 
experimental determination of C’,, it is important to investigate the effect of terms 
of order B® and B® on the decay constant. Hitherto, experimenters have fitted 
their 2 versus B? curve with a formula of the type (23), since neither the 
theoretical estimate of the higher order terms was known nor the experimental 
accuracy was sufficient. In deriving expressions for terms of order B® and B®, 
we shall put K=O in (15) 7.e.we consider a case with no absorption, since it 
simplifies the algebra and the results of interest are unaffected. The expression 
for the smaller of the two roots of 4 can be put in the form 


4 4 B? 4 4B 
A= |—= ——— Ff = F 
5 V3 pit 9Vx F Vx ve 9Va | 


2 4 2 my) -214 
«fr-[e2 s oa (lon. ~ 18} h-2 * | ; * (28) 


Yu 3yi1 V1 
‘I 71, 
where F= (3 ae is). (26) 


It is now straightforward to show that terms of order B® and B® in 2 are, 
respectively 


2 BS 7. gla iad ee 1 
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In order to estimate the magnitude of these correction terms we shall assume 
hat 1,, is constant. In that case, it is easy to show that 


n=, \Kig=t and i, 41 and F= —3, (29) 


where 1 stands for the transport mean free path and corresponds to A, of 
ca paper I. Equations (27) and (28) then, respectively, reduce to 


- Tek |- “i * a 

‘a 9Vaynil 64] ” ee 

iv 4B8 ] 1 

— and Se (- =) . (31) 
3Vayiz 128 \ 72 


a In the same approximation, the term of order B* from (17) is 


= ay oe 
Viv 8 


(32) 


Adding (30), (31) and (32) and using (19) we have for the diffusion cooling 
- constant the expression 


2. 4 72 SS 
o=(1 ae eee | ie. (33 a) 
84 128x3yi1 16441 % 


The smallest assembly used by deSaussure and Silver in their experiments on 
beryllium had B?=7x10-?cm~*. Remembering that /=1.51 cm for beryllium 
and 4y,,=M,=0.32 em™* (values from paper I), we find that the second and 
third terms in the square brackets in (33a) contribute nearly 16 per cent and 
1 per cent, respectively. It is enough, therefore, to take only terms up to 
order B®. It should be pointed out that the correction again depends sensitively 
on the values of J,,. The reason we have used the constant transport mean free 
path approximation in estimating the magnitude of these corrections is that we 
get a better agreement with the experimental value of C of deSaussure and 
Silver. 

It is instructive to note that the above corrections at low temperatures are 
temperature sensitive, whereas at high temperatures they remain constant for a 
given value of B®. At low temperatures / increases and y,, decreases and, there- 
fore, the correction term increases rapidly in magnitude. To reduce this correc- 
tion the assembly has, therefore, to be very large at low temperatures. 

Since the main correction to C is due to the term of order B®, we shall 
now investigate it in the next higher approximation 7.e., we now retain the first 
three terms in the Laguerre expansion and shall again assume that the transport 
mean free path is constant. We shall give only the final result. Instead of 
(33 a) we now get 


Se ered agente i) ae 
1611 % 8 Vu 4. 192 yoo 
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where we have neglected terms containing y,, since they are very small. On 
comparing (33b) with (33a), it is evident that the correction due to the B® 


2 
term is reduced by the factor (1 ae re pe all (22) ), which is nearly 7/10 if 
4 yon 192 \voe ’ 
we take the value y,.=2y,, (see Appendix). Thus, the 16 per cent correction 


as estimated before will now be reduced to only 11 per cent. Therefore, even 
for the smallest assembly of deSaussure and Silver, the correction due to terms 
of order B® is only 50 per cent of the non-Maxwellian correction. And on the 
average it would not amount to more than five to six per cent. Hence the 
main correction to the diffusion cooling constant is the one which is non-Max- 
wellian. 


5. Estimation of the mean energy of neutrons inside and outside the assembly 
The mean energy per neutron is given by 


fad(r, 2, t)dx 
g=- _9— yo Mit.) 


i okey 58 Np (Tf, t) 
0 


2 


=2-2-a, (34) 
where in deriving (34) we have used (4) and (13). 


For the value of A under consideration, « from equation (14), if we retain 
terms up to order B*, is given by 


B 1 oe | Inn \? 
ga Se E teteee 
il " TA 18 V2 al ee el ; (soe) 


which for constant J; reduces to 


Bl Bod 


= Abe oe = 
6V2y,, 18/2 yu 8 eh 
Using the above value of « in (34) we have 
ae E Bl ee 


Thus for B*=0, @=2, as one would expect. If we now substitute the values 
for B’, | and y,, of the previous section, we find that the last two terms in 
the square brackets in (36) contribute nearly 13 per cent. Therefore, in a finite 
beryllium assembly with B*=7x10-? cm?, the inside mean energy per neutron 
is 13 per cent less than what it would be in an infinite assembly. Again this 
deviation would increase rapidly at low temperatures. ae a 
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i \ : 
where A, it given by Eq. (21). Substituting the value of « from (35a) and 
retaining terms up to B’, we have 


ae Cees 2 BA, (A, 3\ (A, Ai 
4 toa Z| 5 (a) cre al a 
5 
© 2 . A, 
- Thus for B=90, (@)ng= =. 
: Ap 
=2 


if the transport mean free path is taken to be constant. In the case of a be- 

ryllium assembly with B’=7 x10? em?, the correction term in (38) is of the 
order of 10 per cent for constant /,,. Therefore for the smallest assemblies con- 
sidered the mean energy per neutron is nearly 10 per cent less than what it 
would be for an infinite assembly. 
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If we assume that /;, is independent of ener 
from (39b), ®in—Zour=0, %.€, a measurement ae mean 6 y ov 
assembly gives also the inside mean energy. But in ite ty « 

neutron energy and the difference of the energies inside and outside 
given approximately by Eq. (39b). For beryihara this difference is ; 


per cent. However for graphite where aed is almost equal to 2, the difference 


is insignificant. If the moderating a is very small, the term involving 
B? in (39a) should be taken into account. ; 


6. Non-Maxwellian correction ; 


In calculating the non-Maxwellian correction to the diffusion cooling constant, 
we have to retain terms at least up to 7=2 and k=2 in the Laguerre expan- 
sion. We shall neglect absorption since it does not influence the result we are 


interested in and the algebra is considerably simplified. Putting p= (7) and 
x 
remembering that px =p,;, Eq. (7) becomes 
2. a 
p> Dig ae +5 la BS + yy,0° =0, (40) 


where we have made use of Kq. (11). For the case under consideration Eq. (40) 
gives the following three similtaneous equations: 


ds© ds ds? 


Poo aE TPs aa eae, dt ae Es? } 
ee i Be sg + - Leh B* 3s = 0, (41a) 
ds ds ada 


Doses oa dt qe a dt hoy ae dt += 5h Bo 


1 1 
Fg hn Be 8 — Sly, Bes + 1 8 + yo, 8 =0, (41 b) 
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In order that there may exist a non-trivial solution for 6, « and £, the deter- 
m inant formed by their coefficients in (43) must vanish. We would then get a 
cubic equation in 4. Making use of the following: 


4 ee VE? his 


] 
4 Po~ 5 ? Pio 473° Pi Gs | 
(see Eq. 16), and af (44) 
V3.2 js ese esi Vax 
Po 16 ’ Pru 54 16 an. Po T9387 | 


iy 
which relations are easy to prove and remembering that Piz=Pxri we get the 
following equation: 


#+ai?2+bA+c=0, (45) 
where : 
3 _% , 256 E {8 28 pe Pps UL 
, oe 6 ean 3 4 02 2 12 4/2 01 64 00 
q ; 213 
ji 31, Bhat Brut era 97 , (46) 
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1 ' ' 
— Ip; bon) — oe Cs — loo tax} = 9 BY {2 19 (lor tos 


— Too Iya) + 711 (loo t22 — [2) + 7/22 (loo tax = %01)} 


1 7 
+ 3 Blo it ae len v23| : (48) 


To determine 4 we proceed as follows: ~ 

We know that the approximate value of 4 which satisfies (45) is 4,, and it 
is given by Eq. (18) without the first term since in this case we have taken 
2a=9. We write A, as 


4,=A(1—e), (49) 
ES 
where A =—]}., (50) 
3Vx , 
Be EN 
and === (: = TA. ; 51 
311 M00 < 2 V2 Wr 


Let the value of 4, which satisfies Eq. (45) be A,+, where x is small. Then 


A, +a is one of the roots of Eq. (45) in which we are interested. x is then 
given by 


S44 a ft+b 
Rt akitba +e . 
3Ai+2aA,+6 
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i as | “ra (- VB lqg+4ly, +751 he) + Yoo (« Lat slo y2 Sire (55) 
eam (53) terms only up to order B* have been retained. It is shown 
in the appendix that for a heavy monoatomic gas y..=2y,, and y,,=0. For 


an actual solid this will, however, be approximately true. Assuming that y,,.<y,, 
and y;2<72, we get a very simple expression for « and therefore for A. Hence 


2 Loo 2 4 ee 
A= ee (11-5! aa io) 


3Vx OV avn 2Vo 9 Vx Yo 
(: 3, les 
__ 2g B 2 Bt 1 \2 y 02° “g “00 
lng -—=) | 144-5 | 56 
3Vx ~ Ty ( >) : 


The second term in the square brackets in (56) represents the main term of 
the non-Maxwellian correction to the diffusion cooling constant, arising as a 
‘result of considering the first three terms in the Laguerre expansion for the 
neutron-velocity distribution function. It represents an improvement over the 
result deduced from the variational principle in paper I. We shall now estimate 
this correction which we call AC/C for the case of constant transport mean 


- free path. x 
| AC_yy ( (V3/8)?- li, 


2 (= 7) @ 


| 


; i.e., AC/C is nearly 19 per cent. 
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spectrum equilibrium neutron spectrum « 
decribed in terms of a neutron temperature which differs from the moc 
temperature. a i 10a 

It is not difficult, though somewhat tedious, to estimate the correction due 
to terms higher than the third in the Laguerre expansion. We have not done 
so but our guess is that these terms may contribute nearly five to six per cent. 
Thus the total non-Maxwellian correction to C is nearly 25 per cent. The higher 
order corrections can be easily estimated on the basis of a heavy gas model. 


~ nearly 


i “ ° ° o ~ ° 
7. Variation of the diffusion cooling constant with temperature 


The variation of the diffusion cooling constant with the temperature of the 
moderating assembly is of special physical interest in the range of low tem- 
peratures i.c., 7 <6, 6 being the Debye temperature of the solid, since it is in 
this range that the effect of crystalline binding plays a dominant role. The 
region of high temperatures i.e., 7720 is not of much physical interest, though 
it may be of some practical importance. As far as the author is aware, no 
reliable experimental results concerning the variation of the diffusion cooling 
constant with temperature are available. deSaussure and Silver (1958) are 
planning to do these experiments in a beryllium assembly in a fairly wide 
temperature range. Their very preliminary experiment at 100°K indicates that — 
the neutron decay does not become exponential even after a fairly long interval 
of time This obviously means that the time constant of the thermalization pro- 
cess is very large at low temperatures, as expected from the theory. 

It is not possible to derive an analytical expression showing an explicit tem- 
perature dependence of the diffusion cooling constant. We shall, therefore, con- 
sider here the two limiting cases: (i) 0/T’—>co and (ii) 0/70, the formulae 
for which have been derived in the appendix. 


Case (i) 0/T’' oo. 


3 
From Eq. (A. 19) Dye(7) TS, 
T 3 
and from Eq. (A. 13) M,«x (7) , 
i 
And since = % M, we have 
as ib 
ox $ | ee 
x(n) ™ /t*(5) 
9\8 
e o=(n) 7 (57) 


Thus the diffusion cooling constant would increase very rapidly with the de- 
crease of temperature for low temperatures of the moderating assembly. This 
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uld imply that the leakage rate would be so large that probably a thermal 
uilfbrium would never be attained at these temperatures. 


_ Case (ii) 6/T 0, 


Dy= 5 1,6 
=const TJ, 
and from Eq. (A. 12) | M, = const. 
_ Therefore, C= constant. (58) 
1 Equation (58) will probably be true even for 7’~6. Detailed calculations in 


- the case of beryllium indicate that (58) is true in the range 6/T =0 to 0/T =2. 
It is necessary to make some comments here regarding equation (57). Equation 


3 (57) is true only in the limit (7) oo, In deriving (57), it has been assumed 


that almost all the contribution to D, arises from the inelastic scattering otf 
neutrons by the lattice, which would not be true in general for all solids ex- 
cept at very low temperatures. However, in the case of beryllium even at 
- 150°K, the foregoing assumption may be true, since the Brag cutoff for be- 

ryllium is as high as 60°K and as we know, below the Bragg cutoff there is 
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8, 


at these temperatures. Here the contribu D, would st 
probably come from — scattering, which is not temperature sensitive. : 
that case D,«Tt and M,«(7'/6)° as before and, therefore, Cx(6/TP. fy 
"In Fig. 1 the diffusion cooling constant for beryllium is plotted as a function 
of 6/T. These calculations are exact. The non-Maxwellian correction as givet 
by (56), and which is nearly +20 per cent, has not been added. It may 
noted that as (6/7') increases, M, decreases and therefore the correction due 
to terms of order B® as given by (33) also increases and it has an opposite 
sign to that of the non-Maxwellian correction. From the figure it is clear that 
from (0/T)=0 to (0/T)=2, C almost remains.constant; it then slowly increases 
and after (0/7)=4 it starts increasing rapidly. It would be interesting to 
verify these results experimentally. 


2% 


8. Non-diffusion correction 


In Section 3 we have, in the diffusion approximation, derived an expression 
for the decay constant of the fundamental mode in a multi-velocity theory. It 
is, however, well known that even in a one-velocity group theory the decay 
constant would contain terms of order B* and higher, which ordinarily are 
negligible for very large assemblies. As we are mainly concerned here with small © 
assemblies for the determination of C, it would be necessary to estimate the 
magnitude of these terms. We might call this correction the non-diffusion cor-— 
rection, for the simple reason that no diffusion approximation is made and the 
transport equation is solved under certain approximations. We shall neither 
give here any derivation of the results nor the assumptions under which it is 
derived, since it is essentially contained in Davisson’s book (1957) (see also 
Sjostrand, 1959). 

We start with Eq. (5.5) of Davisson: 


1s Lt 


oar Er aa 
where Iie | (3.0) -2) Feae-v / for @av=3,-2; (60) 
0 0 : : 
0-3, /3,-%, (61) 
and ey, (62) 
iL ‘ 
where B* is the geometrical buckling. 
The decay constant 2 is given by 
A Sy +v Dee (63) 
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, it is easy to show that 
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n deducing (64), the scattering is assumed to be isotropic in the laboratory 
rem. The anisotropy in scattering can be taken care of as in Section 17.3 
Davisson’s book. For the moderators of interest like beryllium, beryllium 
oxide and | graphite, it is sufficient to consider terms of order (*) in the dif- 
ferential scattering cross section. In that case one need retain only the P, term 
in the Legendre expansion of the differential scattering cross section. We can 
then replace >, in (64) by Di= >;(1—cos 6), where >;, is the average macro- 
scopic transport cross section. Equation (64) then becomes 


2 B * 
Rese pte (1 tigget )— dat BD, (14 


B 
35n 7 16 Taal -). Ht 


15 
, Since J,, is of the order of unity, it is seen that even for assemblies for which 


_ B?=7x10-? em’, the non-diffusion correction term is nearly one per cent or so. 
From Eq. (18) we see that the multi-velocity correction in the constan trans- 
2. 


= 7 
* port mean free path approximation is eM The ratio of the multi-velocity 
‘a 2 
- to the non-diffusion correction is therefore, 
(AA) multi-velocity 15 1 


(AA) non-diffusion 6 1,,M,’ (66) 


which in the case of beryllium is nearly 5. Thus one can always neglect the 
non-diffusion correction in comparison to the multi-velocity correction to the 


decay constant. 
9. Results and discussion 


Before we proceed to compare the theoretical values of the diffusion cooling 
constant with the experimental values, it is necessary to make a few general 
remarks. For an accurate experimental determination of the diffusion cooling 
constant, the moderator assemblies have to be necessarily small, and in that 
case, according to formula (33), terms of order B® in the expression for the 
decay constant should be taken into account. This correction is negligible only 
if B'l,,<8y,, or 2M,. Hitherto, experimental curves for the decay constant 
have been fitted with the relation of the type 4=3,v+D,B’—CD,B*, thus 
neglecting terms of order B®. mn 4 

In calculating the numerical values of the diffusion cooling constant from 
our formula, we shall assume that the transport mean free path remains con- 
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‘li distribution corresponding to the temperature of the 
ons are smoothed out. Detailed ca ions of Singwi and Kotha 
show that the above assumption is fairly good in the case of graphite 


so good in the case of beryllium and beryllium oxide. With our present kno 


coherent approximation, it is not possible to calculate the quantities Ig, toy, U1, 
ete. with an accuracy of more than ‘five to ten per cent which, however, is 
-not enough for the calculation of quantities like Paes 


he SVB} 
1 a) and (: ——|, ) 
(1, 2/3 02g “oo 


in formula (56). We shall now consider the cases of beryllium, beryllium oxide 
and graphite separately. 


i“ 
‘ 


(a) Berylliwm.—Here we shall consider only the more recent measurements of 
deSaussure and Silver (1958):-Their value of the diffusion cooling constant C 
is 1.1+0.9 cm? and the maximum value of B? in their experiment was 7x 10~” 
cm’. For this value of B’, the correction term of order B® in (33b) is —11 %, 
and the non-Maxwellian correction of formula (56) is nearly +20 %. If, how- 
ever, we assume that on the average the former correction is of the order of 
five per cent, the net correction to the diffusion cooling constant is +15 %. 
Thus the calculated value (C).q is 0.9 em? (with D,=1.25x10* cm? sec?), 
whereas the experimental value, (C).xp is 1.1+0.9 cm®. The agreement between 
the two values is fairly good. It would further improve by five per cent or so 
if we add higher order non-Maxwellian corrections. Whether this agreement be- 
tween theory and experiment. is significant is questionable for two reasons: in 
the first place the experimental error is very large and secondly we have as- 
sumed that the transport mean free path is independent of energy. The experi- 
mental value of C of Antonov et al. (1956) is 2.5-++0.8 em?. In the calculations 
the value of M, as given in paper I has been used. Our value of M, differs 
slightly from Nelkin’s value which is 0.36 cm™', since he has used for the 


Debye temperature for beryllium a value of 930°K whereas we have taken it 
to be 1000°K. 


(b) Beryllium oxide.—So far as this moderator is concerned, the only available 


edge of the vibrational spectra of these solids and with the use of the in- — 


Seiad 


experimental value of C is that of Iyengar et al. (1957). Their value of OD, 7 


is 3.85x10? cm* sec’. The maximum value of B? in their experiment was 
6.0x10°* em*. Using the theoretical values of 1,,=1.58 cm (Singwi and Kothari 
(1958)) and M,=0.20 em~', we find from (33) that terms of order B® contribute 
nearly —12 % for the maximum B?. Assuming on the average this contribution 
to be —6 % and the non-Maxwellian correction to be nearly +20 %, we get 
a net correction of nearly +15 % to the value of © as given by (17) for con- 
stant transport mean free path. Our calculated value of C Dy is 2.0 em! sec“, 
which is nearly half the experimental value. Probably the experimental value 
is too high. The value of D, as measured by Iyengar et al. in the same ex- 
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riment is somewhat lower than the value quoted by other workers (Koechlin 
t al. (1956)). Tyengar et al. have also measured the variation of the diffusion 
cooling constant in the temperature range from room temperature to 140° C, 
and they find that it increases very rapidly with the temperature, at 140° C 
it has increased by more than a factor of three. This result is hard to under- 
stand and reconcile with the theory. In view of all this it would be advisable 
to repeat these experiments before a satisfactory comparison could be made 
_ with the present theory. 
__(¢) Graphite.—Since the first measurements of Antonov ef al. (1956), Beckurts 
_ (1957) has made a careful measurement of the diffusion cooling constant. His 
_ experimental value of C is 7 +1 cm”. The maximum value of B? in his ex- 
_ periments was 5.7x10°? cm~*. Using the value of 2.7 em. for /;, and M,=0.068 
= ecm.* (this value of M, is due Nelkin, 1958), we find that the correction due 
- to terms of order B® is —7 %. In spite of the fact that here (B)ma, is small, 
this correction is comparatively large and this is due to the small value of M, 
for graphite. Taking an average correction of —4% and a non-Maxwellian 
correction of +20%, we find that the calculated value of C is 7.0 em”, in 
very good agreement with the experimental value. Contrary to the case of 
beryllium, this agreement in our view is indeed significant since neither the ex- 
_ perimental error is large nor the assumption of a constant transport mean free 
path in the case of graphite is bad. 
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APPENDIX 


In the incoherent approximation, the formula for the differential scattering 
cross section can be written in the form 


Gn (@ — Wo), (A-1) 


where o, is the free atom cross section, k, and k are, respectively, the wave 
numbers of the incident and scattered neutrons, and fiw, and fw are the cor- 
responding energies in units of hn, Om being the maximum frequency. 2 W 
is the Debye-Waller factor and is given by 


=F (w + @—2 Vere cos 6) 7 (0), (A-2) 
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where @ is the angle between k and kj, M is the mass of the atom and m — 
that of the neutron and y (0) is defined by i 


70) | 
i 


f(€) is the frequency function and Bs ky here is the Boltzmann’s con- 
‘0 


d &; (A-3) 


stant. 

The notation used here is that of Sjélander (1958). The G-functions are de- 
fined in Sjélander’s paper (Eqs. II°22 and II°27). 

For the moderators under consideration, it is a very good approximation to 


retain in (A-l) terms up to the order (F). Then (A-1) becomes 


aoa = oo y/ 0. 
dQdw 42 @o 


[6,4 a) 7 ( Ye Feo) ey (0) (4, — Gy) V wa, cos OF (A-4) 


From (A-4) it is clear that in the Legendre expansion of >, (@)—-@, 2—’) 
only the terms J=0 and /=1 need be retained. This fact was used in deducing 
equation (65) from equation (64) in Section 8. 


Calculation of M,.—In the calculation of M,, we are concerned with only 
inelastic scattering. In the summation in (A-l), n=0 term is, therefore, to 
be omitted. Performing the angle integration in (A-1) and retaining terms 


2 
up to (7) we have 


do  ™m w \* oy (m\? ‘w\? 
ie oY (0) (=) (w@ + @o) G, (@ — wo) — =, (7) y* (0) fl 


(3 w* + 3. wo + 10 wa) {G, (w — wo) — 1G, (@— a )}. (A-5) 


Here, we shall only consider terms of order (3). The expression for M, 
(Paper I, Eq. 4) is 
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Now 
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A-7 
for |@—w,|>1. an 
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— 


oo 6/T 
TV nee 1 ot Se 
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where the lower limit is —e if e<6/T, and 26 if een 


: T 
_ Eg. (A-10) can also be written as 
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1 1 
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8\6 
3 /T\3 co 6/T 
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~ calculate M,, (A—11) has to be evaluated numerically. The only property of 
the solid which enters on the right hand side of (A-11) is the Debye tem- 
_ perature. 
We have used here the incoherent approximation but the difference between 
the value of M, in the case of beryllium as calculated by (A-11) and the 
value obtained by using the one phonon coherent calculation and the two 
_ phonon incoherent calculation is very small. The present accuracy of the ex- 
periments does not warrant a better calculation of M, than that given by 
(A-11). We might add here a remark for graphite. It is well known that a 
_ simple Debye spectrum as given by (A-8) is not at all satisfactory for graphite, 
and neither is the model of Krumhansl and Brooks very satisfactory as shown 


_ by the inelastic scattering cross section cal 
and Singwi (1957). Schofield and Hassitt 


culations for cold neutrons by Kothari 
(1958), however, find that if one uses 
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a Gaussian form for the G@-function in (A-5), one gets a good agreement with 
the measured inelastic cross section. In that case M, can be calculated as be- 
fore except that one uses the Gaussian form of the function G, as given by 


formula IT'41 of Sjélander (1958), where the constant A, can now be adjusted — 


to give the right fit with the experimental values for the cold neutron cross 
section. We have not yet carried out this programme. 

We shall now consider the two limiting cases of (A-11): (i) when the tem- 
perature is high, ie., 7'/@-—co and (ii) when the temperature is low, i-., 
§/T—>co. Case (i) will give us the classical limit for M,. 


Case (i). The main contribution comes from the second term in (A-11). 
We, therefore, have 
M, 3(T sli eh | 
—2-~—|— ta eda, A-12 
‘is ( 5) ee € | x ax ( ) 
0 


since x is small compared to ¢ and e’~1+2. The right hand side of (A-12) 
is equal to unity. Therefore 


M,=4€& do. 


The above result though rigorously true in the high temperature limit, will 
probably hold good for 7'>6. 


Case (ii). 6/T—> oo. 


In this case it is the first term in (A-11) which gives the main contribu- 
tion. We, therefore, have 


M, ei) , i 20? 
= |= te d 
43, i) OO oie 
0 


: 0 


3/a\° 3 
=3(5). (3) raver. (A-13) 
It must be remembered that (A-13) gives the correct temperature dependence 
but as regards the magnitude of M 2, it gives only an order of magnitude an- 
swer. The intermediate case when poh but not very much greater is of 


practical interest, since for the~moderators under consideration this condition 
is satisfied at room temperature and a few hundred degrees above. In this 


ee it is not possible to derive limiting expressions of the type (A-13) or 


Variation of D, with temperature.—By definition 
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"Substituting the above value in (A-14), we have 


B= 5 (33,(2) (F) reac) (2) [ererae /[eterae | 
e a 2 0 i 
=n a 3 
4 =4(82.(#)ra/aca/) -(-2) (7) (A-18) 
a 
_ Hence 
7 9\2 
: Dy« (7.) Ts. (A-19) 
: Evaluation of yy, and y22.—From equation (9) 
: d 3 (L> Le. L© 
Yu= amen da’ {3,(ar2') LP (a) LP (x) « 


= 35 (a! >a) LP (a!) a! e” LP (a)}. 
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Now ’ LY (0) =a? 6246. 
We have for 75, the following expression: ~ 


ay 


Vos = al | dx a {>; (t>2') a et (x? -62+6)? 
— s(x’ >a) a’ en® (uz? —6a+6) (2"? — 62" 4+-6)}. 


Using the condition of detailed balance and remembering ‘that only those terms 


of the integrand which are symmetrical in x and 2’ will give a finite value for 
the integral, we have 


i 
yanzg | [ dxaw >; (a >a’) xe* pee? (e+ x’)? +18 (¢—2')? 
—6(a—2')? e+e), “ (A-21) 
m cre bi 
In the (i) approximation >, (x—2’)-is given by 


3: T 3 , + =e , ; L 
De (era) = 553 (77) (5) (=) = aera | eet feos 


In general (A-21) has to be evaluate 
cerned with the value of the ratio 


d numerically. Since we are mainly con- 
Yo2/11, We shall evaluate (A-21) in the 
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yr w Ee = (A-23) > 


larly one can show that in the high temperature limit y,.=0. In general 
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